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CORRELATIONS BETWEEN SDSS TYPE IA SUPERNOVA RATES AND HOST GALAXY PROPERTIES 
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ABSTRACT 

Studying the correlation of type la supernova rates (SNR) with host galaxy properties is an impor- 
tant step in understanding the exact nature of type la supernovae. We use SNe la from the SDSS-II 
sample, spectroscopically determined masses and star formation rates, and a new maximum likeli- 
hood method, to fit the Scannapieco and Bildsten rate model SNR = AxM + BxSFR, where M 
is galaxy mass and SFR is star formation rate. We find A = 3.5±^ x lO- 14 (SNe/yr)(M )" 1 and 

B = l-3±g;| x 10- 3 (SNe/yr)(M©yr- 1 )- 1 , assuming overall efficiency of 0.5. This is in reasonable 
agreement with other determinations. However we find strong evidence that this model is a poor fit 
to other projections of the data: it fails to correctly predict the distribution of supernovae with host 
mass or SFR. An additional model parameter is required; most likely this parameter is related to host 
galaxy mass. Some implications of this result are discussed. 
Subject headings: supernovae: general, galaxies: fundamental parameters 



1. INTRODUCTION 

Supernovae of type la (SNe la) played a pivotal role 
in the discovery of the accelerating Universe (Riess et 
al. 1998, Perlmutter et al. 1999), and are among the most 
reliable standard candles on cosmological distance scales 
(Sullivan et al. 2011, Barone-Nugent et al. 2012). How- 
ever, a complete understanding of SN la progenitor sys- 
tems has yet to be achieved, even though credible models 
for such progenitors exist (e.g. Hillebrandt & Niemeyer 
2000, Nomoto 1982, Hachisu et al. 1996). Many methods 
have been developed to discriminate between proposed 
models, including the theoretical modelling of progeni- 
tor systems (e.g., Han & Podsiadlowski 2006, Wang et 
al. 2010), and the measurement of delay time distribu- 
tions (e.g., Mennekens et al. 2010, Totani et al. 2009, 
Sand et al. 2012; Maoz et al. 2012; see also Greggio 2005, 
2010). 

One approach used in constraining the nature of SNe la 
progenitors is studying the correlation between SN la 
rates (SNR) and the properties of the galaxies in which 
they form (e.g., Maoz & Mannucci 2011). The environ- 
mental properties studied so far include redshift (e.g., 
Dilday 2010a and 2010b), host galaxy age (e.g., Gupta 
et al. 2011), environment galaxy number density (e.g., 
Cooper et al. 2009), ejecta velocity (e.g., Foley 2012), 
and, most relevant to this work, mass and star forma- 
tion rate of the host (e.g., Sullivan et al. 2006). 

In this paper, we investigate the correlation between 
SN la rates and the mass and star formation rate of the 
host galaxies in which they formed. One popular model 
for this correlation is the " A + B" model (Scannapieco 
& Bildsten 2005), which states that the supernova rate 
is proportional to a linear combination of galaxy mass 
and star formation rate: SNR = AxM + BxSFR, 
where SNR is the SN la rate, M is galaxy mass in 
stars, and SFR is star formation rate. This relation 
leads trivially to a more physical relation SNR/M = 



aoyiQmvic.ca 

School of Astronomy and Space Sciences, Nanjing University 



A + Bx(SFR/M), where SNR/M is the rate per unit 
mass, and SFR/M is specific star formation rate (which 
is a function of star formation history, independent of 
mass). 

Substantial effort has been invested in determining 
the values of A and B since the work of Scannapieco 
and Bildsten (2005). See Fig. [1] for a comparison of 
some of these results. Neill et al. (2006) and Dil- 
day et al. (2008) used volumetric SN la rates to infer 
these parameters. Of particular interest is the work of 
Sullivan et al. (2006), who obtained stellar masses and 
SFRs of each individual host (and field) galaxy by fit- 
ting PEGASE2 spectral templates to multiband pho- 
tometry (Fioc & Rocca-Volmerange 1997, Le Borgne 
and Rocca-Volmerange 2002). They concluded that 
A = 5.3 ± 1.1 x lO- 14 (SNe/yr)(M )-\ B = 3.9 ± 0.7 x 
lO _4 (SNe/yr)(M0yr -1 ) -1 , values that we will return to. 

On the other hand, Smith et al. (2012) used photo- 
metric masses and SFR's to demonstrate that the A + B 
model was not a good match to SDSS-II supernova data 
(Frieman et al. 2008). They proposed an alternate model 
of the form SNR = AxM x + BxSFRv, with x ~ 0.7 
and y ~ 1. This is, on the surface, broadly consistent 
with the mass dependence found by Li et al. (2011) from 
the Lick Observatory Supernova Survey (LOSS) (though 
see §6). 

The primary objective of this work is to obtain fits to 
SNR models using spectroscopically-deterimried masses 
and SFRs from the SDSS DR7 MPA/JHU value-added 
catalog (Kauffmann et al. 2003, Brinchmann et al. 2004, 
Tremonti et al. 2004), making these values somewhat 
independent of, and possibly more reliable than, those 
based on multiband photometric data. This is of sig- 
nificance since it has been pointed out (e.g., Forster et 
al. 2006) that the uncertainty in measurements of star 
formation histories is an important limiting factor in the 
determination of the values of A and B. In addition, we 
use a new fitting algorithm for A and B based on max- 
imum likelihood, which is more reliable than previous 
fitting methods. 
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In §2, we present the data samples we use in our stud- 
ies, which include the sample of SNe la, the sample of 
photometric galaxies, and the subset of spectroscopic 
galaxies (itself a subset of the photometric sample) . Also 
presented in §2 is the essential preprocessing of the data, 
most notably the matching of the SNe la to their host 
galaxies in the photometric sample, and a brief introduc- 
tion to the MPA/JHU masses and SFRs, in addition to 
a comparison with other masses and SFRs. Our results 
for A and B are given in §3, and a detailed examination 
of the fit of the model to the data is given in §4. Finally 
we discuss alternate models in §5, and the implications 
of these alternate models in §6. 

2. DATA 

The Sloan Digital Sky Survey (SDSS, Abazajian et 
al. 2009) Data Release 7 catalog (DR7) contains 357 
million unique objects, of which nearly 930,000 are galax- 
ies for which spectra are available. From these spec- 
tra, masses and star formation rates have been ob- 
tained (Kauffmann et al. 2003, Brinchmann et al. 2004, 
Tremonti et al. 2004, hereafter MPA/JHU); see §2.2 for 
details. Supernova identifications were taken from the 
SDSS-II Supernova SurvejH, which is a 3-year (2005- 
2007) survey conducted within Stripe 82 of the SDSS. 
Stripe 82 covers an area of nearly 300 square degrees in 
a belt along the celestial equator, and contains ^20,000 
of the spectroscopic galaxies mentioned above. For each 
year the survey was conducted, supernova imaging was 
conducted for 3 months, resulting in a total observation 
period of 9 months. 

2.1. SDSS II Type la Supernovae 

From the complete sample of 660 SDSS II supernovae 
(Frieman et al. 2008), we extract 520 spectroscopically 
confirmed SNe la, 503 of which were observed during the 
3 observation seasons in 2005, 2006 and 2007. (For the 
rest of this paper, "SNe" refers to SNe la unless stated 
otherwise.) Our sample is larger than that used by Smith 
et al. (2012), since we include SNe within z = 0.05. We 
omit 17 SNe observed in 2004, since the observation win- 
dows and completeness of these SNe would be hard to 
gauge. 

Dilday et al. (2008, 2010) have shown that the identi- 
fication efficiency t z in the SDSS-II Supernova Survey is 
> 0.7 out to a redshift of 0.15, and >0.5 out to z=0.25 
(e.g. Fig. 8 of Dilday et al. 2010). In fact the me- 
dian z of SN in our (spectroscopic host) sample is much 
lower ((z) = 0.11) than for the SDSS SN la sample as a 
whole ((z) — 0.2) because of the severe magnitude limit 
r = 17.77 for the spectroscopic sample of hosts. Thus 
incompleteness e z — 0.7 appears to be a reasonable as- 
sumption. We also assume an observation window T = 9 
months in length, and detection efficiency e t = 1. The 
size of the SN sample is therefore assumed to be the in- 
trinsic number occuring in the sky within a timeframe of 
€tt z T = 0.5 year. This assumption is discussed further 
in §6.1, but we note here that different values of £te z T 
result in a simple scaling of supernova rates, and hence 
of A and B values. 

2 http://www.sdss.org/dr7/ 

3 http://sdssdp62.fnal.gov/sdsssn/snlist_confirmed_updated.php 



2.2. SDSS Spectroscopic Masses and Star Formation 

Rates 

The SDSS DR7 MPA/JHU value-added catalog con- 
tains derived masses and SFR values for the spectro- 
scopic galaxies, including 19987 galaxies in Stripe 82 
(our spectroscopic sample). The mass estimates follow a 
methodology similar to that of Kauffmann et al. (2003), 
namely a grid search over a library of star formation 
histories (Bruzual & Chariot 2003) for the most prob- 
able mass-to-light ratio. (This is slightly different from 
the approach of Kauffmann et al. (2003), which matches 
D„(4000) and H8a, but it has been shown that the mass 
estimates of the two methods^ do not differ by more than 
~ O.ldex.) The methods used for obtaining SFR mea- 
surements were similar to those proposed by Brinchmann 
et al. (2004). Galaxies were divided into three classes: (1) 
"SF", (2) "low S/N SF", and (3) "AGN, Composite and 
Unclassifiable" . The "SF" class was processed by fitting 
5 emission lines (H a , Hp, O III, N II and S II) to those 
of simulated galaxies in a spectral library obtained using 
the code by Chariot & Longhetti (2001). For the "low 
S/N SF" class, Brinchmann et al. used a simple conver- 
sion factor between attenuation-corrected H a luminos- 
ity and SFR. The "AGN, Composite and Unclassifiable" 
class was processed by means of the empirical relation 
between SFR and D4000. 

One problem with this methodology is that passive 
galaxies are always assigned a (very small) SFR. This 
can be seen in Fig. [2j where MPA/JHU spectroscopic 
galaxies are plotted in the mass-SFR plane. It is demon- 
strated in §4 that this does not significantly affect our 
results. 

To test the reliability of the MPA/JHU masses and 
SFR's, we compare their results with those of the VESPA 
catalogue (Tojeiro et al. 2009), which estimates the star 
formation history of a number of galaxies in SDSS Stripe 
82 using spectroscopic methods. VESPA masses are ob- 
tained by summing the star formation history; VESPA 
SFR's are obtained from a SFH averaged over the red- 
shift range to 0.11. A comparison of these quanti- 
ties with their MPA/JHU counterparts is made in Fig. 
[3l From this figure, it can be seen that the mass and 
SFR estimates from the two catalogues are in reasonable 
agreement with each other, with offsets of 0.25 and 0.24 
for log M and log SFR, and rms scatter of 0.18 dex and 
0.48 dex, respectively. 

At this point it is relevant to ask how much bet- 
ter spectroscopically-determined masses and star forma- 
tion rates are than their photometric counterparts (as 
used by, for example, Sullivan et al. 2006, and Smith et 
al. 2012). The advantage of using spectroscopic data is 
due to several effects. 1. Redshift estimates (required for 
both the SN and field samples - e.g., §3) are much better 
for spectroscopic data. Using the same sample as ours, 
Smith et al. (2012) showed that photometric redshifts re- 
sult in negligible offsets (0.03 dex) but a scatter of > 0.2 
dex in log M and log SFR relative to using spectroscopic 
redshifts. 2. Photometric methods for obtaining M and 
SFR from SDSS data rely on fitting evolving stellar pop- 
ulation models to 5 broad photometric bands. The sim- 

4 http://www.mpa-garching.mpg.de 
/SDSS/DR7/mass_comp.html 
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plicity of these models, coupled with the wide range of 
degeneracies between model parameters, will result in 
larger random errors, particularly for SFR. By compari- 
son spectroscopic methods allow the use of far more data, 
including emission lines, absorption lines, and spectral 
breaks that are better tuned to the determination of M 
and SFR. 3. From a comparison of two different catalogs 
above, we have already seen that there are significant 
systematics in the determination of spectroscopic M and 
SFR. For all the reasons mentioned in point 2, there are 
likely to be even larger systematic errors in photometric 
M and SFR measurements. These systematics are of the 
greatest concern; none of our conclusions are affected by 
random errors in M and SFR. 

2.3. Host Matching 

We adopt the host matching algorithm of Sullivan et 
al. (2006). We use SDSS DR7 r band data to calculate 
a dimensionless parameter i?25 for every potential host - 
SN pair. i? 2 5 is the angular separation of a SN from its 
prospective host, measured in units of the size of the 25 
magnitudes/arcsec 2 isophote, and allowing for the shape 
of this isophote (see Fig. 2} . The host candidate is then 
identified to be the galaxy with the lowest i?25- SNe la 
with redshifts falling more than 3<r from the photometric 
redshift of the host galaxy are then discarded. Using such 
a selection procedure, for hosts with spectroscopy, host z 
and SN z are always in agreement to a precision of 0.001. 

How large can R25 be without introducing significant 
contamination? To answer this question, we conducted 
Monte Carlo simulations with randomly distributed ar- 
tificial SN positions, keeping in mind that a majority of 
contamination is due to SNe for which the host is not 
seen. We find that the contamination is ~ 8% when 
only SN-host candidate pairs with R25 < 3.8 are identi- 
fied as genuine matches. Contamination could be larger 
if we were to take clustering effects into account, but 
such effects are hard to quantify. Applying this crite- 
rion, we find 351 matches for the SNe within the SDSS 
DR7 database. Of these, 53 are galaxies in the spec- 
troscopic sample, with MPA/JHU masses and SFRs. It 
may be argued that since ~ 170 SNe were unmatched, ~ 
8% x 170 ~ 14 SNe of the 351 could be random matches. 
We consider this number to be relatively small in com- 
parison to the 351 matches. Also, this contamination es- 
timate corresponds to the expected contamination from 
using i?25 alone as matching criteria. The additional red- 
shift constraints would help further lower the number of 
random matches. The question of contamination is re- 
examined in §6.2. 

To test the host-matching procedure above, we use 2 
separate methods. The first is to visually examine the 
images of the SN - host pairs, relative to their R25 = 
1 isophotes. The second is to plot the number of SN 
matches versus the R25 criterion adopted (Fig. [5J. We 
can see that there is a flat region beyond R25 = 3.8, 
which we assume to be the regime where random matches 
between SNe and non-host galaxies become dominant. 

We also compare the discrepancies between our host- 
matching method and the "closest angular distance" 
matching method used by many other authors (e.g. 
McGee and Balogh 2010). We find that the two matching 
methods give different results for ~10% of the SNe. 



3. FITTING A AND B 

We adopt a fitting algorithm based on maximum like- 
lihood. Using Poisson probabilities, the likelihoods of a 
galaxy hosting/not hosting a SN la are 

P(S l \SN) = S ie - S % and P(S l \SN) = e" 5 ', (1) 

where Si is the expected SN la rate for each galaxy within 
our observing window, calculated for every galaxy indi- 
vidually by means of the mass and SFR data. (None of 
the hosts in our sample hosted more than one SN.) Mul- 
tiplying these likelihoods for every galaxy, we obtain a 
log likelihood of obtaining our dataset of 

19987 53 

lni = - V] £„ + V] i lnS m) + constant. (2) 

n— 1 m— 1 

where 53 is the number of galaxies hosting SNe, and 
19987 is the total number of galaxies in our Stripe 82 
spectroscopic sample. 

We perform a grid search to find the maxi- 
mum likelihood A and B values, which are A — 
3.5±{$ x lO- 14 (SNe/yr)(M )- 1 and B = 1.3±g;| x 
lO~ 3 (SNe/yr)(M0yr -1 ) -1 . Fig. Uplots the probability 
contours and error bars for our grid search, while Fig. [7] 
compares these results with the observed SNR at differ- 
ent specific SFR (where specific SFR sSFR = SFR/M) 

While exceptionally simple, the maximum likelihood 
analytical method above does have the drawback that 
it neglects statistical uncertainties in the determination 
of M and SFR. The method could be modified to in- 
clude the effects of these uncertainties, however, given 
our conclusion (§4) that the A+B model is not valid, 
we have chosen instead to simply simulate the effects of 
these uncertainties. We have conducted a set of Monte 
Carlo simulations by seeding random, artificial SNe into 
the spectroscopic sample, and have included the effects 
of uncertainties in M and SFR. We find that A and B 
are biased by +0.05 and -0.05 when the effects of uncer- 
tainties in M and SFR arc included. (This bias is in the 
sense "derived value minus true value".) The mass term 
and the SFR term contribute 27 and 26 SNe la respec- 
tively, implying that there are as many prompt SNe la 
as delayed SNe la according to the model. 

Our final results (A = 3.5tg;? x lO- 14 (SNe/yr)(M )" 1 
and B = 1.3±o'j x 10- 3 (SNe/yr)(M Q yr- 1 )- 1 ) are con- 
sistent within error bars with most of the literature. For 
a summary /comparison of these results, see Fig. [TJ 

4. INCONSISTENCIES IN THE A + B MODEL 

The A and B values above are optimal in a maximum 
likelihood sense. However, there is no guarantee that the 
A + B model provides a good fit to our data; more tests 
are needed. 

As a test of the A + B model's consistency, we rank the 
spectroscopic galaxy sample of 19987 galaxies by mass, 
calculate the cumulative percentage distribution of SNR 
for both the A + B model and the observations, and 
apply a Kolmogorov-Smirnov (K-S) test to the two dis- 
tributions (Fig. [8]). The A + B model with our derived 
values for A and B is rejected at the 99% confidence level. 
This result was also found by Smith et al. (2012). 

To investigate the possibility that other A and B val- 
ues could have passed the test, we plot (Fig. [8]) the cu- 
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mulative distribution functions of our galaxies for both 
an A-only model (SNR=AxM) and a B-only model 
(SNR=£> xSFR), keeping in mind that any set of values 
of A and B must fall between these two distributions. 
(Note that for the A-only and B-only models, the values 
taken for A and B do not affect the result of the K-S 
test.) The A-only model is rejected at an even higher 
degree, while the B-only model passes the test (~ 50% 
rejection). This shows that our data do not support a 
larger value of A/B. 

Passing the mass-ranked K-S test is a necessary but 
insufficient condition for a SNR model to be plausible. 
Almost all models pass the sSFR-ranked K-S tests, in- 
cluding the A + B model (Fig. [9]), but the same cannot 
be said for the SFR-ranked K-S tests, for which both the 
generic A + B model, the £?-only model, and all mod- 
els in between are rejected (Fig. fTOj) . This means that 
a smaller value of A/B will not pass the test. Since 
we already know that a larger A/B value will not work 
(see above), we conclude from these tests that no A + B 
model can match our observed data. The results of these 
KS tests are shown in Table [2] 

The KS test is usually valid only when the measure- 
ment uncertainties in the observed sample are small. 
Therefore the (possibly significant) uncertainties in the 
mass and SFR measurements brings the validity of our 
above methods into question. In order to check the va- 
lidity of our K-S test under such conditions, we use a 
Monte Carlo approach. We generate many artificial sam- 
ples of supernovae using the best fit A and B values with 
the sample of 19987 observed M and SFR values, add 
the effects of uncertainties in M and SFR, and generate 
cumulative distributions of log M and log SFR. The re- 
sults are completely consistent with the K-S tests: the 
A+B model is rejected at >99% probability for log M 
and >97% probability for log SFR. 

The maximum likelihood method for fitting A and B 
returns a probability which can be used as a measure of 
goodness of fit. This probability can be interpreted using 
Monte Carlo simulations. Allowing for the fact that this 
probability depends on the number of supernovae gener- 
ated (which varies from run to run), we find that more 
than 99.9% of all simulations possess a peak probabil- 
ity greater than that of the best fit to the observations. 
Thus we conclude, independently of the K-S tests, that 
the A+B model is not a good match to the data. 

To present an alternate, and somewhat more illustra- 
tive, view of this discrepancy, we separate the galaxies 
by mass into two equally-sized groups divided at log 
M=10.7, and plot specific SNR vs. specific SFR (as in 
Fig. 6 of Sullivan et al. 2006). The resulting plot (Fig. 
Hip shows that low-mass galaxies have a systematically 
higher specific SNR than their high-mass counterparts, 
by a factor of 3-5. Applying different mass cuts consis- 
tently results in the same trend. 

A plot of the differential distribution functions of SNR 
predicted by the A + B model vs the actual observed 
SNR (Fig. [P2|) further illustrates the issue, where the 
curves have been scaled to show the relative absolute 
SNR obtained by each model. From Fig. [12l it can 
be seen that both the predictions of our best-fit A + B 
model (green line) and that of Sullivan et al. (2006) (blue 
line) underpredict the rates of supernovae hosted by low 
mass, high specific star formation rate galaxies, while 



overpredicting the rates in high mass galaxies. The red 
line is a (somewhat unphysical) modified A + B model 
with a constant background supernova rate added (§4.1). 
The two (unmodified) A + B models are in reasonable 
agreement with each other. 

To show how large an effect our uncertainties in the 
determination of the A and B parameters may have on 
Fig. [T2l we plot Fig. [T3l which shows the differential 
distributions for the la upper and lower limits for A and 
B. It can be seen that the uncertainties cannot account 
for the discrepancy at the low mass end. 

As mentioned in §2, the MPA/JHU catalog system- 
atically overestimates the SFRs of passive (non-star- 
forming) galaxies, assigning them a small but significant 
SFR. We redo the A + B fits by setting the SFRs of all 
galaxies with an sSFR smaller than 10~ 115 (M /yr) /M Q 
to zero; our results are unaffected. 

5. MODIFICATIONS TO THE A + B MODEL 

We have demonstrated that the A + B model does not 
reproduce our data. To account for the discrepancy, most 
notably the observed SNR excess in low mass galaxies, 
we now try modifications to the A + B model. Table 
[2] summarizes the results for some of the many models 
tried. Due to our small sample size, and also to intrinsic 
uncertainties in the MPA/JHU measurements of masses 
and star formation rates, there exist strong correlations 
between model parameters when more than 2 parameters 
are used. 

As previously found by Sullivan et al. (2006), Ax M- 
and B x SFR-orAy models do not match the observed 
SNR distributions. The AxM + BxSFR + C and 
(AxM + BxSFR)(l + CxM~ 1 ) models are particularly 
interesting, since they both pass the KS tests, and might 
have a physical explanation (§6). 

We also investigate the model proposed by Smith et 
al. (2012): SNR = AxM x + BxSFR v . We attempt to 
obtain best-fit values for the parameters of our own us- 
ing maximum likelihood, but find that our sample is too 
small to overcome the correlations between 4 parameters. 
Using Smith et al.'s values of x = 0.72, and y = 1.01, 
we refit A and B using maximum likelihood; the result 
is A = b.StH x 10- n ,B = l-0±g;| x 10~ 3 (appropri- 
ate units in the Mq - yr system). Using either these 
values, or the A and B values given by Smith et al. 
(A = 0.41 x lO- 10 ,^ = 0.65 x 10" 3 ), we find that the 
model passes the mass-ranked test, but fails the SFR- 
ranked test at 99% confidence. A very small value of 
x = 0.06 ± 0.10 fits the data better, and passes both KS 
tests - this is probably a manifestation of the C term 
mentioned above. 

6. DISCUSSION 

In this paper, we have obtained A and B values from 
spectroscopic data. However, further analysis reveals 
that there is an intrinsic inconsistency between the data 
and the A + B model: mass-ranked KS tests reject all A 
and B values with a higher A/B ratio than our best-fit 
values, while SFR-ranked KS tests reject those with a 
smaller A/B ratio; furthermore, from Figs. 151 [TT1 and [T2l 
we can see that there exists a SNR excess in low mass 
galaxies. §6.1 discusses some possible issues with the 
data and its interpretation; §6.2 interprets our results. 
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6.1. Data Issues 

We have not applied any completeness corrections to 
the supernovae. The completeness of the sample has been 
discussed in §2; incompleteness will affect the numerical 
values of A and B by less than 0.1 dex, but will not affect 
A/B, nor any of the conclusions regarding the A + B 
model (unless there is strong redshift evolution in the 
SNR - something that is not observed - e.g., Perrett 
et al. 2012). To check this, we redo the work with a 
redshift cut of 0.20 instead of 0.25, and find no significant 
differences. 

Is it possible that the detection efficiency is different 
in galaxies with high star formation, and concomitant 
higher extinction? We have already noted that the de- 
tection efficiency of our sample is higher than for SDSS 
SNe as a whole (§2.1). In fact high specific SFR implies 
lower M and luminosity (Brinchmann et al. 2004) and as 
a result the star forming galaxies (log sSFR > —10.5) in 
our sample have an even lower mean z (about 0.07); this 
means that incompleteness is even less of a problem. A 
simulation of an extreme efficiency variation with sSFR 
shows that, while B will be biased, the main conclusion 
regarding the inapplicability of the A + B model remains 
unaltered. 

The precision of the SDSS MPA/JHU masses and star 
formation rates may affect our conclusions. We have al- 
ready mentioned the effects of statistical uncertainties in 
the derived M and SFR quantities. We repeat our work 
using VESPA masses and star formation rates. The re- 
sulting values of A and B conform better to those found 
by Sullivan et al. (2006) (see Table [T]), but we are still 
able to reject all A + B models using KS tests. See Table 
[3] for a summary of our results obtained using VESPA 
data, and Fig. [T] for our best-fit values for A and B 
obtained using VESPA. 

Our simplified treatment of the observation windows 
(§2) could affect the derived values of A and B (and C, 
where relevant) by a common factor (which we estimate 
to be < 0.1 dex). Again, this does not affect conclusions 
regarding the rejection of A + B models. 

6.2. Physical Implications 

The (somewhat unphysical) SNR = AxM + 
B x SFR + C model is a good fit to the data, but implies 
a constant rate per galaxy independent of galaxy prop- 
erties. Could the progenitors of SNe arising from the C 
term (~ 68±g 7 % of the total number of SNe, la errors 
quoted) be due to a diffuse intergalactic stellar popula- 
tion, not accounted for in MPA/JHU masses? Sand et 
al. (2011) find that about 17% of the SNe la in rich 
clusters occur in an intergalactic stellar population, a 
number that rises to almost 50% in groups (McGee and 
Balogh 2010). As discussed in Sand et al., these num- 
bers are in general accord with the direct intergalactic 
light measurements of Gonzalez et al. (2005, 2007). If 
some supernovae belong to an intergalactic population, 
then nearby galaxies (with approximately the same red- 
shifts as the supernovae) could be incorrectly identified 
as hosts. A detailed modelling of this effect is beyond 
the scope of this paper; however, a C term that requires 
>50% of SNe la to be intergalactic seems implausible. 
An alternate interpetation of the constant C model can 
be obtained simply by writing the rate in the more phys- 



ical form SNR/M = A + B{SFR/M) + C/M. In this 
case C/M provides some mass modulation of the A + B 
model, which we have already noted. 

We now turn to this mass dependence of the A + B 
model, as seen directly in Fig. [11] and also in the fits 
to the two M -dependent models in Table [2] An effect 
similar to this has been seen before in the Lick Observa- 
tory Supernova Survey (LOSS - Li et al. 2011): SNe la 
from LOSS possess a specific rate SNR/M oc M -0 - 5 
(rather than the expected SNR/M oc constant). In the 
case of LOSS, however, reliable SFRs were not available 
due to a lack of spectroscopic data, rendering it difficult 
to disentangle the mass-dependence of mean galaxy age 
or specific star formation rate (e.g. Schiminovich et al. 
2007), which could produce some or all of this effect. In 
fact, a simple model incorporating ages from Gallazzi et 
al. (2005) explains the LOSS data rather well (e.g. §3 
of Kistler et al. 2011), because low mass galaxies tend 
to be younger, and have higher specific SFR's, than high 
mass galaxies. 

Mass-dependence of the simple A + B model has also 
been inferred by Smith et al. (2012), who combine SDSS 
SN la data with photometric masses and SFR's to derive 
SNR/M oc AxM-°- 3 + B(SFR/M). As already noted, 
however, this model fails the SFR-ranked KS test when 
their parameters are assumed. 

In this paper the mass-dependence of the A + B model 
is demonstrated in a manner that is different from, but 
complementary to, the methodology of Smith et al.. 
First, we use spectroscopically-determined masses and 
star formation rates, which presumably should be more 
accurate (though resulting in a smaller sample size). 
Second, we visualize the mass-dependence in a non- 
parametric manner from a plot of SNR/M vs. SFR/M 
for two different mass samples (Fig. [TT]) . Note that 
this is quite different from the SNR/M vs M plot in 
Li et al. (2011), which includes galaxies with a wide 
range of SFR's, and hence may be affected by the mass- 
dependence of star formation rate (though see below). 
Third, we use the fact that an acceptable model must 
successfully reproduce the observed distribution of su- 
pernovae with mass and SFR, a fact that rules out most 
models that do not depend on mass. 

How is this mass dependence of the A + B model to 
be interpreted? The obvious culprit is metallicity, due to 
the well-established correlation between metallicity and 
mass for galaxies (Tremonti et al. 2004, Gallazzi et al. 
2005 and references therein). Based on the environmen- 
tal dependence of SDSS SNR, it has been suggested by 
Cooper et al. (2009) that the SNR of galaxies may in- 
crease with decreasing gas-phase metallicity. Kistler et 
al. (2011) note that for a given initial mass, the end- 
point white dwarf is more massive for lower metallicity, 
a fact that may on its own (and ignoring many other fac- 
tors) lead to an enhanced SNR in low metallicity galaxies. 
We note that neither of these results constitutes a direct 
proof of an inverse metallicity-SNR relation. Theoretical 
models for single degenerate SNe la in fact suggest the 
opposite rate-metallicity relation (e.g. Kobayashi et al. 
2000, Langer et al. 2000). 

An alternate, and somewhat more straightforward, in- 
terpretation for the mass dependence relates to the fact 
that the rate of SN la from a burst of star formation 
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(the so-called delay time distribution or DTD) declines 
with time as roughly 1/t (Totani et al. 2008, Maoz et al. 
2011). Referring to Fig. [TIJ the median log mass for our 
low (high) mass sample is 1.7 x 10 10 M© (1.1 x 10 n M©). 
Thus, from the mean age-mass relation of Gallazzi et al 
(2005), and DTD oc 1/t, we would expect the rate to be 
about ~ 1.5— 2 x higher in the lower mass sample (assum- 
ing that all stars within a galaxy form at the same epoch) 
- marginally consistent with what is observed. However, 
we note (Pritchet et al. 2008) that a wide range of evo- 
lutionary models and composite stellar populations (and 
hence mean stellar ages) converge to a narrow range of 
loci in the SNR/M vs. SFR/M diagram in Fig. [TCI 

While the initial mass function (IMF) is generally as- 
sumed to be universal, recent evidence - both dynami- 
cal (e.g. Cappellari et al. 2012) and spectrophotometric 
(e.g. Spiniello et al. 2012, Ferreras et al. 2012) - sug- 
gests that massive galaxies possess steeper IMF's, with 
a deficit of high mass stars relative to low mass stars 
(see also Kroupa et al. 2011). At a deeper level, this 
may well be a metallicity effect. Since SNe la are ex- 
ploding WD's, and since such WD's originate from main 
sequence stars as massive as 8 M Q , it follows that low 
mass galaxies are more efficient (per unit galaxy mass) 
in producing SNe la. A detailed calculation of this effect 
is not yet possible because, among other things, it re- 
quires a knowledge of the efficiency of SN la production 
as a function of stellar mass. However, it is interesting 
to note that, for fixed total mass and a power-law mass 
function dN oc M a dM, varying the power a by ±0.5 (as 
observed) around the Salpeter value of a — —2.35 pro- 
duces a change in the number of SN la progenitors by a 
factor of 2 — 4 x (depending on the mass range used). An 
IMF effect therefore has some potential as a mechanism 
for the galaxy mass dependence of SN la rates. 

Finally, we note the recently discovered host mass ef- 
fect on SN la peak luminosities (Kelly et al. 2010, Sul- 
livan et al. 2010). A somewhat simplistic hypothesis is 
that both this effect, and the host mass modulation of su- 
pernova specific rates that we have found, are connected. 
Both these effects could be caused by the presence of two 
distinct classes of SN la progenitors, whose frequency de- 
pends on host galaxy mass through metallicity and/or 
progenitor mass. 

7. CONCLUSIONS 

In summary, we have studied how Type la supernova 
rates correlate with host galaxy masses and star forma- 
tion rates. Our supernovae and host galaxies were taken 
from the SDSS spectroscopic galaxy sample, which has 
spectroscopic measurements of masses and star forma- 
tion rates. The reliability of the host sample was checked 
by comparing with VESPA data. Our primary results are 
as follows: 

• We find that maximum likelihood values A = 
3.5±o 9 7 x 10- 14 (SNe/yr)(M Q )- 1 and B = 1.3±g;| x 
10 _3 (SNe/yr)(AfQyr" 1 )~ 1 provide an optimal fit 
to our data for a generic A + B model. This result 
is largely consistent with the literature (Fig. [IJ. 

• The A + B model fails to predict the mass- 
dependence of the specific SNR - specific SFR re- 
lation in our sample. We also show that no set 



of A and B can account for this discrepancy be- 
tween the model and the observed data. This result 
holds regardless of any uncertainties in the observ- 
ing windows and completeness. 

• Modifications to improve the A + B model are 
tested. Mass-dependent models perform well in the 
K-S tests. 

A number of explanations of these results exist. Most 
promising is the galaxy mass dependence of mean stel- 
lar ages, possibly coupled with variation in the slope of 
the IMF with galaxy mass. Further work is needed to 
quantify the effect of these two mechanisms on super- 
nova rates. 
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Figure 1. Results of our fits to the A + B model with error bars 
(filled circle) plotted against other values (squares) in the literature 
(Scannapieco & Bildsten 2005, Neill et al. 2006, Dilday et al. 2008, 
Sullivan et al. 2006, Smith et al. 2012). The open circle denotes our 
results obtained using VESPA data, as mentioned in the discussion. 
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Figure 3. A comparison between MPA/JHU and VESPA entries 
for mass (top) and SFR (bottom). Note the constant offset (dashed 
line), as predicted in Tojeiro et al. 2009, which is the result of a 
calibration offset between the VESPA and MPA databases. The 
scatters (standard deviation) for the masses and SFRs are 0.18 dex 
and 0.48 dex respectively. 



Figure 2. Distribution of spectroscopic galaxies in the MPA/JHU 
catalog (grey dots), with the 53 hosted SNe (filled red circles) on 
a logAf - logSFR plane. Also shown for comparison are the error 
bars on the masses and star formation rates of 5 randomly picked 
galaxies. 
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Figure 4. Diagram illustrating the concept of i?25- The ellipse 
is the 25 mag/arcsec 2 isophote for the potential host galaxy, O is 
the centre, the SN la is at B, and A is the point of intersection 
between the ellipse and line OB. i?25=OB/OA. 
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Figure 5. R25 criteria used (x-axis) plotted against number of 
SNe matched (y-axis). R25 = 3.8, the criterion we use, is shown 
by the vertical line. Note that the curve flattens out very rapidly 
beyond this line. 
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Figure 6. Probability contours for the values of A and B on an 
A-B plane. The red, blue and green contours represent la, 2a, 3<r 
probability respectively, while the circle denotes our final maximum 
likelihood value. The error bars show the la and 2a limits of the 
marginalised distributions of A and B. 
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Figure 7. Specific supernova rate as a function of specific star 
formation rate. The black points are for our observed sample. 
Plotted for comparison are the predictions of our best-fit A + B 
model (black line), the predictions of Sullivan et al. 2006 (dashed 
line), and their data (squares). 
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Figure 8. Cumulative distribution functions of mass-ranked K- 
S tests for our best-fit A + B model (green), the A X M-term-only 
model (red), and the B X SFR-term-only model (blue). The green 
and red lines are rejected by the K-S test, showing that our best-fit 
A + B model and all A + B models with a higher A/B ratio do not 
agree with our data. 
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Figure 9. Cumulative distribution functions of the specific-star- 
formation-ranked K-S tests for the various models. The only 
models which were rejected were the mass-only model (red), and 
the SFR-only model (blue). The remaining models (AxM + 
BxSFR,AxM + BxSFR + CxM~ 1 ,AxM + BxSFR + C, and 
(AxM + BxSFR)(l + CxAf- 1 )) all fall within the hashed blue 
area. As described in the text, this set of K-S tests is largely in- 
conclusive. 
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Figure 10. Cumulative distribution functions of SFR-rankcd K-S 
tests for our best-fit A + B model (green), the AM-term-only model 
(red), and the BSFR-term-only model (blue). The green and blue 
lines are rejected by the K-S test, showing that our best-fit A + B 
model and all A + B models with a lower A/B ratio do not agree 
with our data. This, in conjunction with Fig. [8] rules out all A + B 
models. 
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Figure 11. Specific supernova rate as a function of specific star 
formation rate. The open points correspond to our low-mass sam- 
ple, and the filled points our high-mass sample. The vertical error 
bars indicate Poisson uncertainties for our SNe numbers in each 
bin, and the horizontal error bars correspond to bin size. The 
trend that the open points tend to lie above the filled ones exists 
for all mass discrimination criteria we try. Also plotted are the 
predictions of our best-fit A + B model. Our results have been 
scaled to account for observing window issues. 
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log[M/M0] 

Figure 12. Differential distribution functions of supernova rate as 
a function of host mass. The black line corresponds to our observed 
data, presented as number of SNe observed in each respective mass 
bin. The green and blue lines are the predictions of our best-fit 
A + B model and the results of Sullivan et al. (2006) respectively. 
The red line is the prediction of our best-fit AxM + BxSFR + C 
"constant background" model, and the cyan line corresponds to the 
(AxM+BxSFR)(l + CxM~ 1 ) model. The predictions have been 
scaled by the number of bins used within the plotted domain, so 
that the lines are comparable with those of the observed data. Our 
results have been scaled to account for observing window issues. 
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Figure 13. Differential distribution functions of supernova rate as 
a function of host mass. The black line corresponds to our observed 
data, presented as number of SNe observed in each respective mass 
bin. The red and blue lines correspond to the upper and lower la 
limits due to the uncertainties of our results for the values of A and 
B respectively, plotted by setting both A and B to their upper and 
lower limits respectively. The green line corresponds to our best- 
fit values. The predictions have been scaled by the number of bins 
used within the plotted domain, so that the lines are comparable 
with those of the observed data. Our results have been scaled to 
account for observing window issues. 
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Table 1 Sample of 53 Hosts for A + B model fits 



Host ID logM* logSFR* 



587730847691047799 


10.654 


0.834 






587730848501203452 


9.427 


-0.833 






587731 1 73306008275 


11.118 


-0.688 






587731 18511 4350067 


11.011 


-0.663 






587731 185668849933 


11.034 


-0.760 






587731185669505187 


10.557 


0.575 






587731 1 861 95366060 


10.019 


-1.259 




'/J 


587731186197332148 


9.936 


0.241 




3 


587731511537959123 


9.590 


0.147 




1 


58773151 1546806433 


9.850 


0.220 






587731513142214757 


10.489 


-1.020 






587731513142542420 


10.414 


-0.017 






587731513143328955 


10.845 


-0.898 




CO 

03 


587731 51 31 46671 21 5 


11.395 


-0.681 




587731513679675512 


10.226 


0.345 




3 


587731513679610009 


10.534 


0.576 




1 


587731513693569205 


11.017 


-0.566 






587731514222116993 


11.444 


0.326 






587731514227818648 


10.177 


0.222 


CO 


S; 


587731514231816228 


11.104 


-0.774 




3 


587731514232209584 


11.175 


0.471 


T3 
O 


1 


587731172231872692 


11.178 


-0.930 


^> t 




587731172233183337 


10.795 


0.809 


+3 
(3 




587731172767368214 


10.664 


0.919 


CD 


CO 


587731174914786124 


10.569 


-0.229 


t£ 


'3 


587731185121951943 


10.669 


1.115 


l—\ 


3 


587731185126539408 


9.637 


0.269 


i 


CD 
4^ 


587731185129554046 


11.168 


0.085 


O 

CO 


O .5 

u 


587731 1 851 32568670 

c..'L' i \ ■ > i 1 I..'"..' 1 -Ji-;'vK HI i y > 


10.878 


-1.001 


CD 

K 





587731 187278872773 


11.386 


0.940 


u 

a. 


587731512071028897 


10.892 


-1.102 


c 


(ap 


587731512621465722 


10.421 


0.426 


o 


587731513427624076 


9.357 


-0.470 


o 

cn V 




587734305949483196 


11.551 


-0.368 




7 


588015507661783172 


9.366 


-0.655 




588015507672137965 


10.947 


0.782 




u 


588015507677642829 


11.042 


0.752 




. >> 

05 
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CO 




588015508211368141 


10.247 
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E 


CO 
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B 
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& 
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co 
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u 

CD 
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3 
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Ph 
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588015509292319354 


9.103 
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+^ 




588015509293760750 


10.689 


0.644 


CO 
CD 

m 
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-0.946 




588015509811626061 


10.926 
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* M is in M Q , SFR is 


in M Q /yr 
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Table 3 

Rejection Rates of Different Models by VESPA Data 



Model 


(1 - a M ) 


(1 - asFR.) 


AxM + BxSFR 


>99% 


77% 


AxM + BxSFR + CxM- 1 


>99% 


62% 


AxM + BxSFR + C 


39% 


42% 


AxM only 


>99% 


>99% 


{AxM + BxSFR)(l + CxM- 1 ) 


34% 


77% 


BxSFR only 


34% 


>99% 
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